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Abstract . Perhaps the deepest mystery of our accelerating Universe in expansion is 
the existence of a tiny and rigid cosmological constant, A. Its size is many orders of 
magnitude below the expected one in the standard model (SM) of particle physics. 
This is a very welcome fact, namely if we care at all about our own existence and fate. 
However, we do not have a minimally satisfactory explanation for our good fortune and 
for the failing of the SM at that crucial point. To start with, an expanding Universe is 
not expected to have a static vacuum energy density. We should rather observe a mildly 
dynamical behavior 6A{t) ~ 7? ~ H‘^{t) with the expansion rate H. At the same time, 
it is natural to think that the huge value of the primeval vacuum energy (presumably 
connected to some grand unified theory) was responsible for the initial inflationary 
phase. In the traditional inflaton models such phase is inserted by hand in the early 
epoch of the cosmic evolution, and it is assumed to match the concordance ACDM 
regime during the radiation epoch. Here, instead, we consider a class of dynamical 
vacuum models which incorporate into a single vacuum structure A(77) the rapid stage 
of inflation, followed by the radiation and cold matter epochs, until achieving our dark 
energy Universe. The early behavior of such “running vacuum model” (ACDM) bares 
resemblance with Starobinsky’s inflation in the early Universe and is very close to 
the concordance model for the entire post-inflationary history. Most remarkably, the 
inflationary period in the ACDM terminates with “graceful exit” and the large entropy 
problem can be solved. The model is compatible with the latest cosmological data on 
Hubble expansion and structure formation, and at the same time presents distinctive 
observational features that can be tested in the near future. 
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1 Introduction 


When Einstein first introduced the cosmological term, A, in the gravitational field equations, it 
was assumed to be constant and positive [T]. Its value was naturally fixed to be of the order of the 
critical density. The reasons why Einstein introduced that term are well-known0, and at present 
are of mere historical interest. However, there is one thing that remains intact: its tiny value and 
the belief (at least within the standard lore) that A is a true constant of Nature. Certainly we 
cannot disagree about its value since it has apparently been measured [311]. It roughly remains of 
the same order of magnitude as in the time of Einstein, namely it is of order of the critical density. 
Later on a nonvanishing A was used to try to cure various astronomical problems, until it became 
clear that its theoretical status was quite precarious since it was realized that it was far too small by 
any standards in particle physic units. This was the origin of the so-called “cosmological constant 
(CC) problem” [21 [5]. At this juncture people felt it was time to desperately find a theoretical 
explanation for A being very small, most likely zero. It was thought that some symmetry would do 
nicely the job. Supersymmetry |6] (SUSY), for instance, has been frequently cited as of the time 
when Wess and Zumino explored this theoretical possibility - more than forty years ago [7] . 

But nowadays we know that the CC term. A, is nonzero and that SUSY (and for that matter 
any other extension of the standard model (SM) of particle physics) is virtually impotent to explain 
the CC problem in any less troublesome way than the SM itself. Both the SM and the MSSM [8] 
(the minimal supersymmetric version of the SM) have essentially the same acute CC problem. This 
is of course very disturbing, specially for the health of our cherished SM of particle physics, a model 
whose theoretical status became lately significantly augmented after the proclaimed discovery of 
the Higgs boson [9], the last missing piece to “crown” the particle physics puzzle. The bare truth is 
that despite the tremendous effort and final triumph of the Higgs boson finding m, the “excessive 
success” of the SM makes it “so simple yet so unnatural” m- Without extra physics the Higgs 
boson (and hence the entire backbone structure of the SM) stays radiatively unprotected from 
the physics of the high scales, namely those extremely large scales associated to the grand unified 
theories (GUT’s), which - we should perhaps recall at this point - are not just there to provide 
aesthetical ideas for the unification of the gauge couplings, but are also badly needed to explain 
the essential facts of inflation and baryogenesis, for example. 

Quite surprisingly, it is usually said that there is no trace of physics beyond the SM. However 
this is not quite true: in fact, for all its particle physic success the SM is in manifest disagreement 
with the most basic cosmological observations. If we do not find at the moment any obvious 
disagreement at the microphysical level (through e.g. 5a deviations, or more, in some particle 
physics observable) does not necessarily mean that the SM is a truly watertight theory, since it 
makes a catastrophic prediction as to the value of the vacuum energy and hence of the cosmological 
constant (see |2] for a more vivid account), let alone its complete inability to explain or to hint 
at the nature of the dark matter. The upshot is that some fundamental aspects of the SM of 
particle physics turn out to be in blunt disagreement with the standard (or “concordance”) model 
of cosmology - the so-called ACDM model. We could say that by finding the Higgs boson (or 
even “a Higgs boson”, of whatever nature) we have certified experimentally the reality of the CC 
problem. The CC problem, therefore, is no longer a formal conundrum! 

It is thus more than plausible to say that there is, or there must be, physics beyond the 
SM, and also beyond the ACDM! What is not obvious at all, at least to us, is where it lies the 
new physics that should rescue the two standard models, one of particle physics and the other of 
cosmology, from wreckage. We believe at least on one thing: wherever it hides and whatever it be 
the new physics, it should provide a sober and sound explanation for the nature of the vacuum 
energy and its connection with the value of the cosmological constant. And we also believe that 

^See e.g. Ref. [2] for a recent review. 
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this explanation should involve a dynamical cosmological term, namely one capable of tracing the 
history of the Universe from its origin till the present days in a single unified framework. 

While the traditional “explanation” is the existence of a nonvanishing and positive cosmologi¬ 
cal constant, A, whose energy density equivalent, p\ = A/SttG, is of order of the critical density, 
this cannot be a truly convincing explanation, as in fact an expanding Universe is not expected 
to have a static vacuum [2] . A smoothly evolving vacuum energy density p\{t) = p/^{H{t)) that 
borrows its time-dependence from the Hubble rate H = H{t) — taken as the natural dynami¬ 
cal variable in the Friedmann-Lemaitre-Robertson-Walker (FLRW) background — is not only a 
qualitatively more plausible and intuitive idea, but is also suggested by fundamental physics, in 
particular by quantum field theory (QFT) in curved space-time [12] . This more formal point of 
view, based on the renormalization group approach, in which A = A.{H) is a running quantity with 
the cosmic e^ansion, has been introduced in [JS] [H] US] [E] ITT] dS] |T9| , and recently emphasized 
in misniiiijB Furthermore, some recent applications of these ideas for a possible description of 
the complete cosmic history have been put forward in the literature |2l| . These vacuum models 
represent a conceptually new approach that goes beyond the first phenomenological approaches on 
time evolving cosmological constant, cf. [251 EH Ea [281129]. For alternative formulations and recent 
developments along these lines, see [solEIlIM] and references therein. In general, the notion of a 
variable vacuum energy has also been entertained in the literature from different and interesting 
points of view, also from the historical and philosophical perspective [MllMlES]. 

In the last years, and well within the tradition of the old Dirac’s ideas, an independent source 
of puzzling news has generated also a lot of interest. Frequent hints that the electromagnetic fine 
structure constant aem and/or the proton mass might be changing with the cosmic time (and locally 
in space) are reported in the literature [SHET] - for reviews, see e.g. [38ll39]- Theoretical models 
already exist in the literature trying to explain such phenomena, see e.g. [10]. It is tantalizing 
to conceive that, if aem can evolve with the cosmic expansion, all of the fundamental “constants” 
should change in time as well, including the gravity coupling and the masses of all the elementary 
particles. Recently, these ideas have been linked to a possible time variation of the vacuum energy 
as well [niiMiiis], and this opens a line of thought well in the context of the dynamical vacuum 
energy in an expanding Universe [2] • Since the Newtonian coupling, Gw, determines the Planck 
mass Mp = Gj^ ' , a possible time variation of it would be tantamount to say that Mp slowly 
evolves with the cosmic expansion [42] . 

In this presentation of the “running ACDM model” (in fact a class of models) we show that the 
idea of a dynamical vacuum in an expanding Universe can be the sought-for touchstone enabling a 
substantially improved account of the entire cosmological history as compared to the concordance 
ACDM with rigid A term. The content is as follows: In sections 2-3 we summarize the theoretical 
framework for dynamical vacuum models. In sections 4-5 we discuss inflation in the context of 
the ACDM model and address some important thermodynamical aspects of it, in particular the 
large entropy problem. Subsequently, in Sect. 6, we compare the early cosmological behavior 
of the running ACDM model with Starobinsky inflation in alternative formulations. The low- 
energy implications of the ACDM class are explored in Sect. 7, where we confront these models 
with the most recent cosmological data from distant supernovae (SNIa), the cosmic microwave 
background (CMB) anisotropies, the baryonic acoustic oscillations (BAOs) and the input from 
structure formation. We show that the notion of a running vacuum can be, in principle, compatible 
with the current observations, and suggest that traces of their dynamics should be testable in the 
near future. In the final section we provide some discussion and our conclusions. 

■^Recent comprehensive analyses successfully testing a wide class of dynamical vacuum models in the light of the 
recent observations are also available, see |22[I23| . 
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2 Dynamical vacuum in an expanding Universe 


The idea that in an expanding Universe the cosmological term A and Newton’s constant Gm 
could be variable with time is not new and in fact it can be viewed as reasonable and even 
natural. Dirac’s ideas in the thirties [44 1 145] on the so-called “large number hypothesis” were 
seminal concerning the possible time evolution of the gravitational “constant” Gj^. They were 
disputed by E. Teller [46] and further qualified by R.H. Dicke [47] . It also triggered subsequent 
speculations by G. Gamow |48] on the possible variation of the fine structure constant. Since then 
the subject has been in continuous evolution, and more and more sophisticated experiments are 
being designed to monitor the possible time (and space) variation of the fundamental constants - 
see e.g. [38l [39] for reviews. 

While in the old days the cosmological term, A, in Einstein’s equations may not have attracted 
a lot of attention regarding to its potential time variability (setting aside occasional episodes 
where some astrophysical observations had suggested this possibility), it is natural to entertain 
this option in earnest when we cope with the full cosmological context. This is especially so 
if we take into account that A defines the energy density parameter ^ = K/{S'kGn)-, which 
is interpreted as the vacuum energy density of the expanding Universeu. It should perhaps be 
surprising if an accelerating Universe were to carry a static vacuum energy density throughout 
the entire cosmic history [2]. A more natural possibility, which is perfectly compatible with the 
Gosmological Principle, is that A = A(t) and hence p\ = p\{t). 

Erom that point of view, it is instructive to consider the possible modifications that may undergo 
the basic conservation laws if one makes allowance for the time variability of the fundamental 
gravitational parameters Gm and A. The Bianchi identity satisfied by the Einstein tensor on the 
l.h.s. of Einstein’s equations reads V^G^y = 0, where G^j^y = R^j^y — {l/2)g^yR. It follows that the 

covariant derivative of the r.h.s. of Einstein’s equations must be zero as well: = 0, 

where T^y = T^y + g^j^y p\ is the full energy-momentum tensor of the cosmic fluid composed 
of matter and vacuum. Using the explicit form of the Eriedmann-Lemaitre-Robertson-Walker 
(ELRW) metric, the generalized conservation law emerging from this dynamical framework reads 

^ [GNiPm + Pa)] + ^G]sr H (1 + ujm)pm 3 Gn H (1 + uj\)p\ = 0, (1) 

where Um = PmIPm is the equation of state (EoS) for matter and uj\ = pa/PA is the EoS for the 
vacuum. The Hubble rate H dynamics ensues directly from Einstein’s equations. We generalize 
them for the case of dynamical vacuum (and restrict to the flat ELRW metric): 

= SnGipm + pa) (2) 

2H + 3H^ = —87rG{uJmPm + ^APa) ■ (3) 

The overdot denotes derivative with respect to cosmic time t. One can easily check that only two 
of the equations ([I|)-(l3|) are independent. Eor example, substituting (l2|) in ([T|) we arrive at 

H + in Gn (1 -|- uJm)Pm + dvr Gn (1 -|- 0J\)pA = 0 . (4) 

One can check that this equation can also be obtained by combining (l2|) and ([3]). 

®Ten years later after A was introduced by Einstein [T] to insure a static non-evolving Universe, G. Lemaitre [4^ 
introduced a nonvanishing A to discuss his dynamical models of the expansion of the Universe, strongly motivated 
by E. Hubble’s observations prior to their publication in 1929 m- A few years later, in 1934, Lemaitre discussed 
for the first time [51] the interpretation of the CC term as vacuum energy and its associated negative pressure - see 
e.g. [541EU for further historical discussions. 
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Let us note that the previous formulae remain valid if we sum over all cosmic components 
(matter and vacuum). However, in the particular but rather common situation where there is 
a dominant matter component (e.g. cold matter or relativistic matter), it is possible to obtain 
the evolution law for the Hubble function solely in terms of the vacuum term and that dominant 
material fluid. A simple calculation from the above equations leads to 

3 

H + — (1 + LUm) = 4:7T G {iOjn — coa) pa ■ (5) 

As advertized, in this equation we cannot sum over the matter components, as uim in it stands for 
the EoS of the dominant one. This equation will be useful in the next sections. 

Up to this point we have assumed that the EoS of the vacuum component is completely gen¬ 
eral. While the previous equations account for the general situation (applicable even if ua would 
correspond to a general DE fluid), we shall henceforth adopt the simplest scenario corresponding 
to the vacuum state, i.e. wa = —1, as being a characteristic feature of the EoS of vacuum even for 
time-evolving pA = PA{t)- 

In view of the above considerations, we next consider the following cosmological scenarios 
beyond the concordance model: 

Scenario I: pA = PA{t) is assumed variable, and Gn =const. In this case, Eq. ([T]) implies 

Pm H 3(1 -|- U]rn)HPm = PA ■ ( 6 ) 

Since we now have pA 7 ^ 0 it means we permit some energy exchange between matter and vacuum, 
e.g. through vacuum decay into matter, or vice versa. Obviously if pA = 0 we recover pm + ^ H ( 1 -|- 
'^m)Pm = 0, i.e. the standard covariant matter conservation law. Its solution in terms of the scale 
factor is well-known: 

Pm{a) = (7) 

Scenario II: pA = PA{t) is again variable, but Gn = G^it) is also variable. In contrast to 
the previous case, here we assume matter conservation in the standard form ([7]). As a result the 
following conservation law ensues: 


{Pm + Pa)Gn + GnPa — 0 . (8) 

In this setting the evolution of the vacuum energy density is possible at the expense of a running 
gravitational coupling: G / 0. 

Scenario III: Here we keep pA =const., but Gn = GAfit) is again variable. Now we find: 

GAf{Pm + Pa) + GAf[Pm + 3H(1 -|- Um)pm] = 0 . (9) 

In this case matter is again non-conserved and the gravitational coupling is running. Despite the 
vacuum energy is constant in this scenario, such situation can mimic a form of dynamical dark 
energy since it implies a different expansion rate |53] . 

The above three generalized cosmological scenarios differ from the concordance ACDM model, 
but can stay sufficiently close to it if we consider the recent history of our Universe. Let us finally 
note that the above dynamical vacuum models can be extended at high energies for a successful 
explanation of the inflationary Universe |24] . 

3 ACDM: a vacuum model for the complete cosmic history 

Thus far we have sketched some feasible frameworks for the time evolution of the the vacuum energy 
and the gravitational coupling. Let us however note that with the generalized conservation law ([T]) 
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and Friedman’s equation ([2D is not possible to determine the solution of the cosmological equations. 
We need some information on the evolution laws p\ = pK{i) or G = G{t). Such information has 
been provided in the past based on some phenomenological ansatz (see e.g. [2ZII2SI Ea EB] and 
references therein). However, here we wish to focus on a class of models motivated by the theoretical 
framework of QFT in curved spacetime, see [laiiaillKITlIIHlIISllSOlET] -cf. also the review [2] 
and the long list of references therein. 

Recall that in particle physics we have theories such as QED or QCD where the corresponding 
gauge coupling constants gi run with an energy scale p, i.e. gi = gi{p). The scale p is usually 
associated to the typical energy of the process. Following the same line of thought we can think 
of PA as an effective coupling sensitive to the quantum effects and thereby running with an energy 
scale p representative of the cosmological evolution. In the context of QFT in a curved back¬ 
ground, where gravity is a classical external field, the quantum effects that shift the value of pA 
are exclusively produced by the loops of matter fields m- Naturally p should be fixed on physical 
grounds, but it is quite reasonable to assume that the running of pA should be associated with the 
change of the spacetime curvature because if there is no change in the background geometry there 
is no cosmological evolution at all. In the FLRW metric, this means that p should be related with 
H and its time derivatives. We denote this association in the simplified way p ^ H (recall that H 
has dimension of energy in natural units), but we will be more precise soon. 

In view of the foregoing, the following form has been proposed for the renormalization group 
(RG) equation for the vacuum energy density of the expanding Universe [2[ I14j: 


dpAjp) 

dlnp2 


1 

(dvr)^ 


BiMf + X] 

i i i 


+ ... 


( 10 ) 


In this expression, M* are the masses of the particles contributing in the loops, and Bi,Ci,.. are 
dimensionless parameters. The RG equation (HOI) provides the rate of change of the quantum 
effects on the GG as a function of the scale p. Provided we are interested only on the dynamics 
of the current Universe, we may cut off the series at the quadratic contributions, i.e. only the 
“soft-decoupling” terms of the form ~ Mf p^ will be of significance. Notice that the Mf terms 
are absent, as they would trigger a too fast running of the GG term. As a matter of fact these 
effects are ruled out by the RG formulation itself since only the fields satisfying p > Mi are to 
be included as active degrees of freedom. Being p = 0{H) (as indicated above) it is obvious that 
such condition cannot be currently satisfied by the SM particles - see, however m- The leading 
effects on the running of p\ are, according to Eq. (llOj) . of order Mfp^ ~ Mf and hence it is 
dominated by the heaviest fields. In the context of a typical GUT near the Planck scale, these are 
the fields with masses Mj ~ Mx ^ Mp. For instance, in Ref. m a specific scenario is described 
which is connected to the effective action of QFT in curved spacetime. 

Let us recall that because of the general covariance of the effective action, among the list of 
possible terms emerging from the quantum effects one expects only terms carrying an even number 
of time derivatives of the scale factor. If expressed in terms of the Hubble rate (which is the most 
convenient quantity to parameterize the extra contributions), it amounts to terms of the form 
H, H^, H^H etc. In contrast, the linear terms in H (and in general any term with an odd 
number of derivatives of the scale factor, such as H H, H etc) are not expected since they 
would be incompatible with the general covariance of the effective action [laiiHiiisiiiniiii]. Let 
US remark that at low energies only the and H terms are relevant. The higher order ones can 
however be important for the early Universe [MIIMIISSIISSIET]. 

As we have agreed, p ^ H is the natural association of the RG-scale in cosmology. However, 
a more general option is to associate p'^ to a linear combination of and H (both terms being 
dimensionally homogeneous). Adopting this setting and integrating (IlOp up to the terms of 0{p'^) 
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it is easy to see that we can express the result as follows: 


Pa{H, H) = ao + aiH + a2H‘^ + 03 H , (11) 


where the coefficients Oj have different dimensionalities in natural units. Specifically, oq has di¬ 
mension 4 since this is the dimension of pAj and 02 have dimension 2; and, finally, 03 , 04 and 
05 are dimensionless. The “running vacuum Universe” (ACDM) is the extension of the ACDM 
model based on a dynamical vacuum energy density of the form (lllh , stemming from the basic RG 
equation (jlOl) . While higher order term are still possible, that expression contains the basic terms 
up to four derivatives of the scale factor, and hence encodes the basic potential of the model both 
for the low and the high energy Universe. 

Let us now consider a particularly simple and illustrative case of the running ACDM Universe. 
Suppose that rather than associating p? with a linear combination of and H we would just 
set p? = (in this case the linear combination reduces to just one term and we can just adopt 
the canonical choice p = H.). In this situation we have oi = 03 = 05 = 0 in (llip . The remaining 
coefficients can be related immediately to those in ( 1101 ) . and one can show that the final result can 
be cast as follows: 

^ + if) ’ 

where cq has dimension 2 and we have introduced the dimensionless coefficient u and the dimen¬ 
sionful one Hj. Comparing with (fTUI) it is easy to see that 



Es< 


Mf 

Ml ■ 


(13) 


The dimensionful coefficient Hj absorbs any other dimensionless factor, but it is unnecessary to 
further specify its structure since it represents a physical quantity (connected with the mechanism 
of inflation) and can be determined (or at least bounded) by observations, as we shall see in a 
moment. But let us first discuss the interpretation of the dimensionless coefficient ly. The sum in 
(|13h involves both fermions and bosons. Coefficient v plays the role of the /3-function coefficient 
within the structure of the effective action in QFT in curved spacetime. This is confirmed by 
the fact that depends on the ratio squared of the masses of the matter particles to the Planck 
mass, which is indeed the expected result in particular realizations of the RG in curved spacetime 
- see e.g. |18] . As we shall see, H[ stands (to within very good approximation) for the Hubble 
parameter in the inflationary epoch, which is the only epoch where the higher order term 
can be of relevance. During the inflationary period H =const. {H = 0), so at least in the pure 
inflationary stage the terms we have dropped from Eq. (|lll) should not be determinant. These 
terms, however, can be important for the different modalities of reheating, just after inflation. 

Obviously the dynamical vacuum model (I12p aims at providing an unified description of the 
entire cosmic history, valid from inflation to the present days. We will confirm if this is the case in 
the subsequent sections. For the current Universe it is enough to consider (fT^ up to the quadratic 
term, or type-Al vacuum model: 

PAiH) = (co + (Type Al). (14) 

Note that if and Hq are the current values of the vacuum energy density and the Hubble 
parameter, then 


Svr G]\[ 


,0 
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Taking into account that at low energy the terms ~ H can be equally important as the ~ ones, 
one may also consider the slightly extended model variant 

Pa{H) = (co + i>H + (Type A2), (16) 

where P is a new dimensionless coefficient; it stems from the original (dimensionful) oi in Eq. (1101) 
after a convenient redefinition. Although we use the same symbol cq, it is understood that the 
connection of cq in ()16p with the current cosmological parameters is not exactly the same as in 
m but it can be obtained easily. 

We will also introduce briefly two more vacuum types (B1 and B2) for comparison. These are 
defined as follows: 


Pa{H) = 


(Type Bl) 

(17) 

Pa{H) = 


(Type B2), 

(18) 


where e is a new dimensionless coefficient and Hq is the current value of the Hubble parameter. 
The dynamical vacuum models could provide an alternative explanation for the dynamical DE in 
the Universe within the context of QFT in curved spacetime, thereby representing an alternative 
to quintessence and other exotic DE options. In the case of B1 and B2, however, the presence 
of the linear term ~ H (rather than ~ H) is more phenomenological because we do not expect 
linear terms in H in the effective action for the reasons explained above. Still these terms can 
be admitted (they could represent bulk viscosity effects [58], for instance) provided cq 7 ^ 0 and/or 
the standard terms (with an even number of derivatives of the scale factor) are also present. In 
the case cq = 0 the above models do not have a well-defined ACDM limit (i.e. none of them 
has in this case a behavior near that of the concordance model), and one can show that this 
is problematic [59| . Recent analyses have confronted these models to observations [221 [23] - see 
also [ISl dT] [ 6 OI EU ESI E3] for previous studies on a variety of possible scenarios. In Sect. |7| we 
provide a summarized presentation of the main results for the most recent analyses. 


4 Inflation in the ACDM model and Grand Unified Theories 

Let us now concentrate on the high energy implications of the unified dynamical vacuum model. 
To simplify our discussion we will adopt the canonical form (|12p since it already contains the 
main features^. The Hubble function for the dynamical vacuum models under consideration can 
be derived from solving the basic equations of Sect. 2. However, these equations depend on the 
particular scenario (I,H and HI) we choose. For the present study we will focus on Scenario I of 
Sect. 2 and hence on Eq. ([6]). 


4.1 Solving the model in the early Universe 


Assuming that there is a dominant matter component in the early Universe (typically radiation, 
= 1/3), we may combine equations (|5]) and (|12l) and we arrive at the following differential 


equation for H: 




l-v- 


cq 



(19) 


^The analysis of the modifications introduced at high energies by the more general structure HU will be presented 
elsewhere. 
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Now, if Eq. (H is to describe an approximate CC term in the late Universe, the constant cq must 
be of the order of the current value of the CC, or Aq — 3co to be more precise, and at the same 
time the term ~ u can only represent a small (dynamical) departure from it, thus |z^| <C 1. 
This is not only what we expect, but in fact what we find when we compare the low-energy models 
(H and (I16p with the precision observational data collected from the cosmological observations, 
we find indeed |z^|, |f'| < 10“^ (cf. Sect. 7). 

In the early Universe, when H is large, we can dismiss the term cq/H'^ <C 1 in Eq. (|19p to a 
very good approximation. This implies the existence of a constant solution in that epoch, namely 

H = Vl^Hi c:^ Hi , (20) 


which corresponds to a de Sitter phase driven by a huge value of the Hubble parameter Hj. We will 
corroborate that Hj is the characteristic value of the Hubble function of the high energy phase of 
the early Universe, i.e. the effective value of H in the inflationary period, from our analysis of the 
energy densities below. In what follows we will neglect u in all practical considerations concerning 
the early Universe. Notice that when H ~ Hj, the ~ term in Ea. (|12l) is suppressed with 
respect to the ~ by precisely a factor |i/| <C 1. In the current Universe, where H = Hq ~ 10““^^ 
GeV (the present value in natural units), the ~ H‘^ term in Eo. ffT^ is suppressed with respect to 
the ~ one by more than 100 orders of magnitude! ~ cf. Eq. (j24h below. 

The value of the inflationary scale Hi can be estimated from the anisotropies of the cosmic 
microwave background as follows. According to the CMB observations, the (primordial) spectrum 
V(^ of the curvature perturbation [M] at typical “pivot” length scales (of a few tens to hundreds 
of Mpc) is ~ 5x 10“®. On the other hand, the theoretical calculation of Vq is usually performed 
in the context of the inflaton model with effective potential I4ff, and yields 


n(^o) 


8 (Ueff) 

e 


( 21 ) 


where (Uefi) is the vacuum expectation value of the potential and e is the standard slow roll 
parameter |64j . Recall that e is related to the tensor-to-scalar fluctuation ratio r = ut/ rig through 
r = 16e. Furthermore, the value of the potential during inflation is related to the GUT scale Mx 
through (Uefr) Mjr. For definiteness and simplicity we take them equal. Obviously this also 
defines the critical density at the GUT epoch, pi = M^, which is related to Hi by means of 


PI 


2,H‘j 
Stt Gn 


a4. 


( 22 ) 


From the previous equations we immediately find: 


Mx = (Ueff)'/" = ' MpVl^^ 



2.2 X 10^® GeV, 


(23) 


and 


Hi 


I Sir M‘jr 

V 3 Mp 


/ r \i/2 

va2/ 


1.1 X 10^^ GeV. 


(24) 


Notice that we have normalized the tensor-to-scalar ratio to the approximate value r ~ 0.2 fur¬ 
nished by the BIGEP2 collaboration [65] - with the understanding that we are waiting for a future 
confirmation/reanalysis of this important experimental result. Notwithstanding, it is pretty clear 
that owing to the fourth and squared roots involved of the factor r/0.2 in the previous equations, 
the final estimates on Mx and Hi are not critically sensitive to the precise value of r. 
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It is remarkable that the obtained value for Mx in (1231) lies in the ballpark of ~ 10^® GeV, 
which is the expected order of magnitude in all viable GUT’s compatible with the current limits 
on proton decay. At the same time, the value of obtained from (1241) is seen to satisfy 


IIl 

Mp 


< 10-^, 


(25) 


which is the condition that the fluctuations from the tensor modes do not induce GMB temperature 
anisotropies larger than the observed ones. Indeed, the vacuum fluctuations become classical a 
few Hubble times after horizon exit, and have a spectrum determined by the Gibbons- 

Hawing temperature Tgh = Hj/2Tr. To be more precise, the primordial tensor perturbation hij 
has the spectrum 


Vhiko) 


QAtt [HiV 


16 

vr \Mp) 


(26) 


Therefore, by requiring that < 5 x 10“® we retrieve the bound ()25p within the correct order 
of magnitude from Vh alone. 

Of course the previous considerations are just an estimate inasmuch as they are partially 
based on the standard inflaton picture, which is not exactly the one we are subscribing here. They 
nevertheless lead to a reasonable estimate, for we could have just reversed the sense of the argument 
and started from Eq. (I26p rather than from (j2ip . Following this logic, we could have next applied 
the bound (1251) - imposed by the measured anisotropies of the GMB - and derived an upper limit 
on the inflationary value of the Hubble parameter directly from (I26h , which is under the constraint 
10 and then immediately derived the corresponding GUT scale (1231) . The advantage 
of this way is that we know that the general form of the tensor spectrum, Pf^ {H/Mp)"^, Is only 
sensitive to the absolute value of the Hubble parameter squared - measured in Planck units - at 
the inflation phase, a fairly general fact that holds good Irrespective of the underlying details of 
the inflation dynamics. In this sense, the previous considerations are essentially safe in so far as 
Eq. (I25p is compatible with them. 

Far away from the inflationary period (i.e. H <C Hj) we find another de Sitter solution of 
Eg. (fTOl) . namely H = [co/(l — whereby A « 3co ~ Aq (recall that \v\ <C 1). This Is 

of course the solution we have mentioned before, which leads to the late time (approximate) 
cosmological constant behavior, i.e. the current DE epoch. Let us emphasize that the pure de 
Sitter phase of this late time solution is achieved only in the remote future, but the fact that v is 
small but not exactly zero Is responsible for a slow vacuum dynamics ~ - cf. e.g. Eq. (|14p - 

in our recent past (still going on at present) which could mimic a dynamical DE scenario. This 
is precisely the situation when v ^ Q matters, and it can affect the present observations. We will 
address the details later (see Sect. 7), but first we further analyze the inflationary phase. 


4.2 Achieving “graceful exit” in the ACDM model 


Let us now elucidate the cosmology of the ACDM model in the early de Sitter epoch. We can 
virtually set cq,v —>■ 0 and assume w = 1/3 since the subsequent matter epoch is presumably 
relativistic. It is convenient to trade the cosmic time for the scale factor through d/dt = aHd/da. 
Eq. (I19p takes on the form 

+ (^1 - = 0 , (27) 

where prime denotes differentiation with respect to the scale factor. The solution of the previous 
equation renders 


H{a) 


Hi 

[l + Da4]i/2 • 


(28) 
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Substituting the above result in Eq. 
in terms of the scale factor: 


m we get the explicit form of the vacuum energy density 


PA{a) = PI 


1 

[l + Da^f ■ 


(29) 


The matter energy density then ensues after inserting the two previous results in Friedmann’s 
Eq. ([2]). We obtain: 


Pr{a) = Pi 


Da^ 

[1 + Da^f ■ 


(30) 


In the previous formulas pi is the same quantity defined in (1221) . We conhrm now that it represents 
the primeval energy density for a —)■ 0, i.e. in the inflationary period, as it is obvious from (I29p . 
Furthermore, since the (relativistic) matter energy density ^ 0 as a —)• 0, it follows that 
ultimately pi is “pure vacuum energy” stored in the very early Universe. As shown by the above 
equations, in the course of the fast evolution the vacuum energy transforms into matter, so the 
Universe becomes eventually dominated by the material component. The maximum density of 
radiation achieved is p^^^ = p//4. Beyond that point the radiation density starts to decrease and 
at the same time the inflationary period (which is powered only by the vacuum energy) comes 
naturally to an end. This is tantamount to saying that the ACDM model incorporates natural 
inflation with “graceful exit”, as we shall further discuss below. 

In the initial period, <C 1, the solution (I28|) can be approximated by the constant value 
(j20p . i.e. H ~ Hi, the vacuum energy density remains essentially constant and the scale factor 
increases exponentially. This fact can be verified from ()28p by integrating from some initial (un¬ 
specified) scale factor up to a value a{t) well within the inflationary period, in which D 1 

still holds good. The result is just the inflationary solution 


a{t) = Qi exp {Hit} . 


(31) 


The meaning of t in this equation is the elapsed time within the inflationary period, and a* is the 
purported value of the scale factor at the time when inflation started (defined to be at t = 0), 
hence a(0) = a*. What happened before that instant of time, we don’t care. We will probably 
never know since the pre-inflationary Universe is not accessible. 

Once the inflationary phase is left behind, the term D can be comparable to 1 or much 
greater. Here the integration of Eq. (1281) must be done with the full expression, and we find 

[ —\l + Dd^Y^^ = Hit. (32) 

Ja, a 

In this case t is (at variance with the previous case) the time elapsed after approximately the end of 
the inflationary period, indicated by t*, and we have defined o* = a(t*). The integration constant 
D is fixed from the condition H{a^) = H*, and it satisfies Da^ = (H//H*)^ — 1. 

We can now show that the “graceful exit” from the inflationary phase can be successfully 
accommodated in this class of models [23]. In fact, considering the limit S> 1 of equation (1321) 
we find 

a ~ , (33) 

which confirms our contention that the vacuum phase decays into a radiation-dominated Universe. 
We can reach the same conclusion from the analysis of the energy densities in the transition 
period. In effect, from (j29p and (|30p we see that in the limit D S> 1 the decay law of matter is 
the characteristic one of radiation density, pr ~ a“^, whereas the vacuum energy density decays 
much faster, pA ~ , so it becomes suppressed as pa/P r ~ once the radiation regime has been 

settled. Such suppression is welcome, of course, as the success of BBN (Big Bang Nucleosynthesis) 
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Figure ll The transition from the inflationary period of the running vacuum model ca into the FLRW radiation epoch. 
The vacuum and matter energy densities are normalized with respect to primeval density pj , and the scale factor is normalized 
with respect to the transition one a* (see the text). The figure shows an initial (and constant) dominance of the vacuum phase, 
and its subsequent decay into relativistic particles (whose energy density rapidly increases) until the Universe is dominated 
by radiation and the vacuum energy becomes negligible. As a result inflation ends with “gracefully exits” into the standard 
radiation phase. We have used SttG = 1 units. 


is not jeopardized in this model. Recall also that when y 0, we have |z/| <C 1 (cf. Sect. 7) so the 
is also harmless at the BBN epoch. 

As a concrete illustration, the evolution of the primeval matter and vacuum energy densities 
from the very early Universe to the present day is highlighted in Fig.[Tl where we can clearly see 
the transition from the de Sitter vacuum-dominated phase into the radiation dominated phase. 


5 The horizon and entropy problems in the ACDM 


From the analysis of the previous section it follows that the ACDM Universe (|12l) starts with a 
huge vacuum energy density of the order ~ - confer Eqs. (I22p and (I23h ~ and therefore without 

an initial singularity, and has no horizon problem. It follows that a light pulse beginning in the 
remote past at t = ti > ti (i.e. shortly after inflation started at U) will have traveled until the end 
of inflation, tf, the physical distance 


dnitf) = a{tf) f 
Ju 


dt' 


a{t') 




aiHi 


(34) 


where a{t) is given by Eq. (13111 during the inflationary stage, and H = H[ remains constant under 
inflation. In practice H need not be strictly constant, it is enough that it remains approximately so 
for some interval ti < ti < tf. The integral (I34p diverges for ti —>• ti if the initial time of inflation 
is a sufficiently early one (viz. if Oj —>■ 0). This proves the non-existence of particle horizons in the 
running vacuum Universe. Equivalently, and perhaps more transparently, the above integral can 
be written as 


r°-i 

dniaf) = a{tf) / 

J ai 


T da' 




Hi \ai ajJ \ai J ^ 


(35) 
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Being ai at time ti > t, exponentially smaller than aj at the end of inflation, we have afjai ;$!> 1 
and so the above integral (hence the horizon) can be as big as desired. The result (1351) is of course 
the same as (fMl) because ai = Oj . 

Notice that x ~ d-n/o^f — /o-i is the angle subtending the horizon from the point where 
inflation starts. This angle is a comoving coordinate, and if it is fixed appropriately at time ti 
(i.e. if it is made large enough such that it subtends at least the current horizon) it will remain 
so for the entire cosmic history^. The exponential solution triggered by a period of approximately 
constant vacuum energy (starting from an early point ai) makes this arrangement possible and 
hence there is effectively no horizon in the ACDM model. 

In contrast, in the concordance ACDM model the above integral gives (e.g. in the radiation 
dominated epoch) du (this is again du ^ H but now H is not constant) and as a result 
—>■ 0 for a —)• 0, which is tantamount to say that the observers become isolated in the past 
(namely, their horizon becomes infinitely smaller than the size of the Universe for t —)• 0 and hence 
cannot exchange information with it) . This is the well-known horizon problem of the concordance 
model, which as we have seen can be overcome in the ACDM model. On more physical grounds, 
we can say that in the ACDM model the local interactions (the exchange of information) are able 
to homogenize the whole Universe since the very beginning. 

5.1 The large cosmological entropy 

We can reconfirm the absence of the horizon problem in the running vacuum model from a ther¬ 
modynamical perspective, namely from the point of view of the so-called “entropy problem” in 
cosmology. Succinctly formulated reads as follows: how did the universe manage to start with an 
extremely low entropy value and develope later a huge one that gave rise to the arrow of time in 
accordance with the second law of thermodynamics? There is no solution to it within the ACDM. 

Let us start by recalling some basic thermodynamical facts of our Universe that are intimately 
connected with the horizon problem. The current horizon is defined by Hq^ and we normalize 
the scale factor as oq = 1- We can take as a fiducial volume of the observable Universe (up 
to an irrelevant numerical factor common to all our estimates). From standard thermodynamical 
formulae of the thermal history |66] we can estimate the total entropy of our Universe multiplying 
the entropy density by its fiducial volume. In natural units, we find: 

So = ^ gsfi T^o - 2.3h-3l087 ~ iqSS . (36) 

Here T^q ~ 2.725°K ~ 2.35 x 10“^^ GeV is the present CMB temperature. Ho = 2.133/i x 10“^^ 
GeV (with h ~ 0.67) is the corresponding value of the Hubble function [1], and the coefficient 
Qsfl = 2-1-6 X (7/8) ~ 3.91 is the entropy factor for the light d.o.f. today, which involves 

the well-known ratio Ty^/T^o = (4/11)^^^ of the current neutrino and photon temperatures. 

However, in the context of the AGDM concordance model, the huge number (j36p cannot be 
understood without generating a phenomenal causality problem, and hence recreating the horizon 
problem again. The reason is that in the ACDM the total entropy contained in the horizon at 
earlier times is much smaller. In fact, we can project the result (j36p backwards in time taking into 
account that it evolves as (T 77“^) . Now, for an adiabatic evolution T ~ a~^. Moreover during 
the nonrelativistic epoch H ~ so that altogether this renders S ~ = (1 -|- We 

conclude that at recombination (zrec — UOO) the total entropy was a factor ~ ~ 10“® smaller, 

or, to be more precise: S^ec ~ 10®^. It means that the current Hubble volume should contain some 
hundred thousand causally disconnected regions at recombination! This is completely unacceptable 

®For this it suffices that a//ai be only ~ 1% of the total expansion ao/af ~ 10^® after inflation. Eqnivalently, 
the number of inflationary e-folds (e^ = af/ai) must be > 26 In 10 ~ 60. 
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since the smoothness of the CMB must have a causal explanation in terms of interactions that 
propagate at subluminal velocities. Such an unsettling situation worsens more and more when 
we travel deeper and deeper into the past, where the number of causally disconnected regions 
keeps increasing inordinately. It goes without saying that this a serious drawback of the standard 
picture, in fact one of the biggest cosmological conundrums of the ACDM model! 


5.2 Solving the entropy problem with running vacuum in a generic GUT 


The situation in the running vacuum Universe is completely different. We shall show that the total 
entropy contained in the current horizon can be computed from Sq Soo Hq ^ , where S^o 
is an asymptotic value attained by the entropy per comoving volume (or “comoving entropy” for 
short) in the early Universe. Specifically we find that S{a) ^ ^ S^o for a S> Ogg, where a^q 

is the vacuum-radiation equality point. This point is placed in the (still incipient) radiation epoch, 
and in it the energy density equals that of the decaying vacuum. We will show that the comoving 
(hence intrinsic, volume-independent) radiation entropy increases very fast {S oc a®) owing to the 
primeval vacuum decay until reaching that saturation value S'oo- It follows that the amount of 
entropy contained in a physical patch of volume V ~ H~^ (taken at any a 3> a^q) increases just 
proportional to the fixed number S'oo times the physical volume of the region, i.e. S ~ Soo V. In 
particular, for the region encompassed by the current horizon, Vq ~ we have Sq ~ Soo . 
Remarkably, this value is of the correct order of magnitude of the entropy today, Eq. (I36|) . Since 
the comoving entropy S/V ~ Soo is conserved after we have left well behind the vacuum-radiation 
equality point, there is no causality problem at all in the running ACDM model. 

To substantiate these claims, let us hrst estimate the equality point a^q- It can be obtained 
from equating the densities (l2^ and (pO]l . i.e. pr{aeq) = p\{aeq), hence Da^q = 1. This determines 
D = a~q and from (l28l) we see that at a = Oeq the Hubble function reads H^q = Hj /\/2 (obviously 
below Hj). As we shall see, the main results of this analysis do not depend on the details of the 
GUT framework that rules the dynamics of the vacuum. However, the precise value of Ueq is of 
course model dependent. We need not know it accurately, but it is convenient to have a numerical 
estimate. Inserting H(,q = Hi/y/2 and a = Oeq in Friedmann’s equation relating the evolution of 

leq up to now, yields 

, 9 *\ ^ -^7U / y* \ -^7U 

^eq - I ) rr^ = ( ) nh ’ (37) 


radiation from the time when H = 


45 J y/MpHi 

where in the last step use has been made of Eq. 


15 y Mx ’ 

This result is of course not exact (as one 
cannot reach the equality point assuming that the evolution is always within the radiation epochb 
but it provides the estimate a^q ~ 3 x 10“^®, which is sufficient as an order of magnitude valuel£| 
for a typical GUT defined at the scale Mx ~ 10^® GeV. 

To compute the entropy generated from the primeval vacuum decay we have to estimate the 
temperature of the radiation heat bath that is formed during the conversion of the vacuum energy 
into relativistic particle^. This post-vacuum heat bath will mimic the reheating process of the 
standard inflationary scenarios, but the process to achieve it is actually quite different. Let us call 
the radiation temperature of the heat bath T^. Equating the radiation density (j30|] emerging from 
vacuum decay in our model to the general form of the radiation energy at an effective temperature 
Tr{a) we obtain: 


Prifl) — Pi 


{a/a, 


eq) 


[1 -b {a/OeqY] 


412 


_ rp4 


a 


(38) 


®The early equality point Oeq is specific of the ACDM Universe, and it should not be confused with the standard 
equality point at which the energy densities of radiation and cold dark matter become equal: ueq — 3 x 10“"^ (i.e. 
at redshift zeq — 3300). The latter occurs much later and is coincident in both models. We have, ueq ■ 

^See [57] and [BS] for recent alternative formulations, and [55] for early developments. 
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Here counts the total number of effectively massless degrees of freedom (d.o.f.) in the heat 
bath at the given temperature. For the standard model (SM) of particle physics, = 106.75, if 
we include the top quark and the Higgs boson since all of them can be relativistic d.o.f. at high 
temperatures. In general gf* will be larger than the SM value when we operate at a GUT scale. 

From Eq. (1381) we find Tr as a function of the scale factor: 


Tr{a) 


Tx 


a/a, 


eq 


[1 + (a/aeg)^] 


411/2 


(39) 


Here 


Tx 


30 pi 


1/4 


30 


Mx 


(40) 


is the value of the temperature associated to the GUT scale Mx through Eq. (|22ll . The temperature 
at a = 0 is T = 0. It then raises (linearly with a in the beginning) until a maximum value, which is 
attained at a = a^q- It is natural to think of such maximum value as a kind of effective “reheating 
temperature” after vacuum decayH: 


Trh 


V2 



Mx- 


(41) 


We point out that, in the conventional inflaton scenario [66] , the reheating temperature of the 
radiation heat bath is not just determined by the initial GUT temperature Tx- A new parameter 
comes into play, the inflaton width and the reheating temperature is Trh ~ g^ ~ 

Tx^JT^h where Hj ~ g, T^^/Ml- 

Notice that neither Tx (nor Trr of course) are very sensitive to the precise value of gf*, even if 
the number of d.o.f. changes by one order of magnitude, because diOli - dlT]) depend on the quartic 
root of < 7 *. Therefore we will just take the SM value g^ = 106.75 as a fiducial estimate for our 
calculations. In this case we find e.g. Tx — 0.41 Mx, which is smaller than Mx but not far 
away from it. As the GUT scale Mx is some three orders of magnitude below the Planck mass, 
Mp ~ 10^® GeV, we are well reassured that our working regime is perfectly compatible with 
a semiclassical description of QET in curved spacetime (which is the precise context where the 
running ACDM cosmological framework is formulated m)- 

At the maximum temperature (|4ip the radiation density is also maximum, so at that point 
the Universe becomes maximally populated of relativistic particles emerging form vacuum decay. 
Eventually both the radiation energy and the temperature decay (for a /§> Ugg), respectively as 

pAa) ~ P, {^y ~ Mi {^y (42) 

and 

Tr{a) ~ Tx — ^ Mx — ■ (43) 

a a 


It follows that, in the asymptotic adiabatic regime, the thermodynamic behavior of the running 
vacuum model boils down to the standard form, namely pr ~ 1/a^ and T ~ 1/a. Thus, the 
conventional BBN picture can proceed normally since the changes affect only the inflationary 
epoch. 

^Strictly speaking there is no “reheating” in this context since the Universe is just heating up progressively from 
a = 0 till the point a = a^q, where it reaches the maximum temperature m- Later on, when the evolution enters 
fully the radiation epoch, the temperature falls down standard as in (113). 
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Having found an expression for the effective temperature of the primeval radiation heat bath, 
the comoving radiation entropy of the relativistic particles that now populate the Universe can be 
estimated. Starting from 


5^ (a) 


Pr{a) +Pr{o) 3 
Tr{a) 


^ Pr{a) 3 
2,Tr{a) “ 


27r^ o o 

—g^T^{a) , 


(44) 


and then using (|39p we arrive at: 


S{r) = ‘^9* {T^alq) fi^) 

We have defined the following function of the ratio r = ajaeq- 


(45) 


(46) 


Since /(O) = 0, /(I) = 2“^/^ ~ 0.35 and /(r) ^ 1 for r 1, it follows that the radiation 
entropy (1451) evolves monotonically very fast (S' oc a® in the beginning) starting from zero, then 
goes through the equality point a = a^q (where is still increasing) and finally reaches the maximum 
value 

2-k^ 

S(a) ^ Soo = ^ 9* {T^alg) (for a > Oeg), (47) 


where it saturates. This substantiates our claim on the rapid increase of the comoving entropy 
until achieving an asymptotic value Soo- 

But of course we still have to prove that we can use Soo to predict the total entropy contained 
in our present horizon, Eq. (I36p . We are now in position to complete this task. First of all let 
us rewrite the expression appearing in parenthesis in the asymptotic formula (|47p in terms 

of the temperature, T^ad, and the scale factor, Urad, deep in the radiation epoch, i.e. a^ad 3> a^q- 
From Eq. (l3^ we find: 


Trad = Tx 


O^rad! 


eq 


[1 T {cLrad/^eq}^] 


411/2 


~ Tx 


^eq 

^rad 


(48) 


It is important to realize that this approximation is valid for any value of arad well after we have 
passed the equality point Ogg, i.e. when we have entered the adiabatic regime of the evolution. 
From (j48P it follows that Tx aeq — Trad ^rad hi very good approximation in this regime. Inserting 
the result in Eq. (|47p , we find that the total entropy in the current physical horizon is predicted to 
be 

27r^ 27r^ 

Soo = — {g, Tl, alj = — (5.,o «o) , (49) 

where in the last step we have used the entropy conservation law of the adiabatic regime, which 
states that 9*T^ad^rad ~ 9sfiT^oCio- Recall that in our normalization oq = 1, and hence we have 
reproduced the result Sq given by ([361) . 

The upshot is significant: it turns out that the theoretical prediction (I49p of the running ACDM 
model for the total entropy enclosed in our current horizon is precisely the observed entropy in 
our Universe, i.e. the huge entropy number (j36p - in natural units. This result may be viewed as 
a potential solution to the cosmological entropy problem in the context of the running ACDM. As 
we have seen, we can assume some generic GUT at the very high scale, which triggers the decay 
of the huge vacuum energy and generates the large entropy in the radiation phase. Remarkably, 
the important CMB constraint (I25p is preserved and the final result for the entropy is universal, 
meaning that it does not depend on the details of the GUT. The universality of the prediction is 
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a reflex of the vacuum decaying dynamics and of the entropy conservation law, which holds good 
in the Universe once the cosmic evolution has entered the adiabatic regime. 

The vacuum model (|12p that we have studied here (in which the highest power of the Hubble 
rate is H^) is just the simplest implementation of the running ACDM as a candidate model for a 
complete description of the cosmic history from the inflationary times to the present day. However, 
the main results are maintained if an arbitrary even power (n > 2) of the Hubble rate is used 
for the higher order term (as required by the covariance of the effective action m)- At the same 
time, one has to keep the additive term so as to reproduce the standard ACDM model at low 
energies. In Sect. [3 we shall see that the observations are compatible with keeping also the H^- 
term of (jl2p . This is interesting as the latter endows the current cosmological vacuum with some 
mild dynamics (as a kind of “smoking gun” of the entire vacuum decaying mechanism) and can 
be interpreted as dynamical dark energy. 

6 Starobinsky inflation versus running vacuum model 

In this section we wish to elaborate on the fact that the class of vacuum models (1121) realize in an 
effective way the Starobinsky type of inflationary regime [70^1 . This occurs somehow through the 
dominance of the highest power in the early Universe, which is of the order of the Starobinsky’s 
correction term B? ~ in the effective action, see Eq. (I50p below. In actual fact the situation is 
a bit more complicated since at the level of equations of motion the term B?' in the action gives 
rise not just to a pure contribution but to a fairly involved combination of powers of H and its 
time derivatives. Still there are some interesting similarities worth noticing. Some of them have 
been recently pinpointed in |54l I56j . and explored in great detail in m- 

6.1 The standard Starobinsky action 

Let us recall that the original Starobinsky model m is based on the following action: 

S = J + bB^^ + S^^tter , (50) 

where b is the dimensionless coefficient of the higher order derivative term B^. It is usually written 
as 6 = Mp/QMg , where Mg is a mass dimension parameter - playing the role of scalaron mass in 
the original model [70]. Recall that B is, if written in the context of the FLRW metric, a linear 
combination of and H, and therefore B? involves terms of the form and H, all of 

which are roughly speaking of order ~ and hence one may envision a kind of close connection 
of the Starobinsky model and the ACDM model. This is true in part, but it is not quite so, at least 
not for the original form (I5flp of the Starobinsky action. In the next section we briefly mention how 
to improve the connection of the running ACDM with an alternative Starobinsky-like action [18] . 
At the moment we want to continue with the action (|50p and show the form of the inflationary 
solution that emerges from it, as this will enable us to compare with the inflationary scenario 
derived from the ACDM model in Sect.4. 

The equations of motion for the case where the energy-momentum tensor of matter contains 
a mixture of various fluids V = 1, 2... with densities and pressures {pn,Pn) can be obtained after 
varying the action (j50p with respect to the metric. The result reads as follows: 

+ 327rG6(V^V,R - g^^^DB + BB^,, - ^B^) 

= STrG'^lpNgiiu - (PN+Pn)UI^U^] . (51) 

N 

®For recent and past detailed studies of Starobinsky inflation, see e.g. [Zll[22], and references therein. 
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Figure 2‘. The inflationary solution of (I55I I corresponding to the Starobinsky model (I50II . On the left it is shown the 
exponential growth a ~ ^ of the scale factor and its stabilization into the radiation regime a ~ We have taken 

Hj = Mp in (I56I I. and b = 10® from the mentioned bound in the text. On the right we depict the corresponding behavior of 
the Hubble function and (in the inner window) the characteristic oscillations when the Universe leaves the inflationary phase 
and enters the radiation epoch in the form a ~ + oscillations. In that window we have set h = 100 to make the oscillations 

more clearly visible. Time has been rescaled as t = (VOOtt/Mp) t, and H = {l/a)da/dt is the correspondingly rescaled Hubble 
function. 


For 6 = 0 we recover, of course, the standard Einstein’s equations. However, for 6 7 ^ 0 the result is 
more complicated. Setting (/i, i^) = (0,0) and = (i,j) and using the flat FLRW metric, we 

find two independent equations: 

^ '^PN + 96TrGb{H^ - 2HH - 6H^H) (52) 

3 


and 

H^ + — + 32TrGb{2H + 12HH + ISH^H + 9H^) = ^Pn- (53) 

If we now just project the result for a single matter component {pr,Pr), and assume that this 
component is relativistic (j)r = Pr/ 3), we can combine the two previous expressions to obtain a 
single differential equation for the Hubble rate: 

2ff^ -hff-h 48TrGb(2ff + Uffff + 24^2^ + = 0 . (54) 


It is pretty clear that this equation (a third order, nonlinear, differential equation) is considerably 
more involved than the hrst order equation (|19p for the ACDM model in the early Universe. This 
shows that the initial analogy between the two models, based on the fact that ~ was an 
exceeding simplification since the models have indeed to be compared at the level of the equations 
of motion. 

Expressing (15411 in terms of the scale factor it becomes a fourth order differential equation: 


a aa 
+ 3— 7 ^—h 


a2 


a a? a? 


- 5- 


,aa^ 


a 


_ p 

3 967r6 


+ 


= 0 . 


This equation is to be solved under initial conditions, which we take as follows: 


a(0) = Oj d(0) = Qi H] d(0) = a* a (0) = ai Hj 


(55) 


(56) 
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and assume Hj ^ Mp in this case. An analytic solution of this problem is not possible in general. 
However, for 6—7-0 the dominant term is the last one of (j55p . and the solution becomes a r\j tl/2 
corresponding to the radiation epoch. In contrast, for large 6 the last term of (I55p becomes 
negligible and, then, the solution is easily seen to be the exponential one: a{t) = The 

numerical solution of Eq. (j55l) under the initial conditions (|56p is given in Fig.[2j Remarkably 
an approximate solution of (1551) . or equivalently of ([MP . can be found to describe the end of 
the inflationary period. Since H remains essentially constant until we are very near the end 
of the inflationary phase (see the straight line in the plot on the right in Fig. [2]), we can solve 
(|54P by neglecting H/H^ <C 1 and all higher derivative terms. In this way we are left with the 
equation STOttGat b H = —1, whose immediate solution is H{t) = Hj — Mpt/576TTb, and hence the 
corresponding scale factor reads: 

a{t) = tti 1152776 . (57) 

Obviously, 6 > 0 in order to have a stable inflationary solution until the inflationary phase termi¬ 
nates at around t/ — 1152 7r6/Mp. The larger is 6 the longer is the inflationary time. One can 
show that if we impose the limit (|25P from the CMB anisotropies the parameter 6 has a lower 
bound of order 10® - see [72]. 

Needless to say the approximate solution (I57P does not perfectly interpolate from the inflation¬ 
ary epoch to the radiation epoch and hence, in contrast to the situation with the ACDM model, 
we do not have an analytically rigorous description of the graceful exit. The solution (j57p is only 
indicative that the inflationary process comes to an end, and that there is a chance for a correct 
transition to the a ~ radiation epoch. But it misses the nontrivial details accounting for the 
reheating stage. Unfortunately, the issue of graceful exit is still a complicated and unsolved matter 
in Starobinsky inflation, despite it has been discussed in different places in the literature, see e.g. 
m and m in different formulations. A more detailed study, though, is still needed m- 

In point of fact, we see from the comparison of Figures 1 and[2] that the connection of the 
Starobinsky model and the ACDM model is not that close after all, at least using the original 
action (|50p . It turns out that the connection with the ACDM model is much better accomplished 
if one adopts the conformally invariant formulation of the Starobinsky model, which is only broken 
by quantum effects; namely, the so-called anomaly-induced effective action (see next section). The 
reason is that only in this formulation the ACDM model can be motivated from an effective action, 
which is certainly not the one given by (j50p . The latter is, in effect, not conformally invariant. 

6.2 The anomaly-induced effective action 

A generalized form of Starobinsky’s inflation, based on the effective action of anomaly-induced 
inflation was considered long ago in uniiiiiiini - and references therein - but the first explicit 
connection of this framework with the running vacuum model ACDM was signaled in |18j . One 
hnds that by taking the masses of the fields into account, using the conformization procedure of the 
cosmon model |76] , the inflationary process automatically slows down and can therefore favor the 
graceful exit. Let us briefly mention the solution of this alternative form of Starobinsky inflation. 
We follow Ref. [18|, where more details are provided. Rather than starting from the action (|50p 
one takes the classically conformally invariant higher derivative action of the vacuum: 

S'vac = J {aiC'^ + a2E + a^V'^R} , (58) 

where ai, 2,3 are dimensionless coefficients, is the square of the Weyl tensor and E is the Gauss- 
Bonet topological invariant in 4 dimensions [l2| . This action can then be completed with the 
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conformally invariant realization of the Einstein-Hilbert (EH) term: 

r _ 

S%H = J [Rx^ (dx? ] ■ (59) 

Here x is a background scalar field that realizes the conformal symmetry at the high energy 
scale A4, presumably close to the GUT scale Mx- The same field and scale are used to conformize 
the masses of all the matter fields in the action |18] following the prescription of m- After X 
acquires a vacuum expectation value of order Ai, that symmetry is spontaneously broken and we 
recover the standard EH term. Now, of course, the conformal symmetry of the total action is 
also intrinsically broken by quantum effects, not only the traditional ones from the trace anomaly 
associated to the 5vac part [12] but also by the additional contribution form The anomalous 

induced part of the action (i.e. that one breaking conformal symmetry from quantum effects) 
follows from solving the functional differential equation: 




(5f 

T-1- 


1 




6r 




= biC^ + b2E + + 


f 

WttGnM'^ 


[Rx" - eidxf 


(60) 


The one-loop values of the /3-functions 6 i, 62 j ^3 are well established since long time ago [12] and 
depend on the matter content of the model. The new ingredient is the /3-function of the conformal 
EH term, which is given by |18j 


- _ 1 Np m?p 1 Nv My 
' M], 


It involves contributions from fermions (F) and vector bosons (V) since the scalar part vanishes in 
the classical conformal limit. The trace-anomaly equation (I60p can be solved following the standard 
method [77], namely from a conformal transformation of the metric extended with 

the new background field: x = e~'^X- this way one finds the corresponding effective action that 
includes the effects of Spjj- [73l[T8]. The relevant part of it can be summarized as follows: 

SeS = J ~ ® (^x)^] + Sraatter + high, deriv. terms , (62) 

where 

Ml{a) = fa) (63) 

is a “running Planck mass” that is sensitive to the dynamics of the conformal factor a. We may 
next fix the conformal gauge as x = Ml (i.e. x = Mle°^) so as to recover the EH form of the action. 
The scale A4 cancels and the resulting dynamics does not depend on it. Finally, we know that the 
FLRW metric is conformally flat in conformal time rj, so that with the choice = a^(r/) we can 
take the flat metric gfj_y = This renders the new curvature terms trivial, R = 0, but of course 
the action depends in a nontrivial way on the conformal factor a = In 0 ( 77 ). Varying this action 
with respect to it and reverting to the cosmic time t one finds a fourth order differential equation 
for a{t) similar to ([5^ . although not identical. We shall omit cumbersome details here [73l IT8] . 
but we remark that an approximate solution of it can also be found in the form (j57|] as follows: 

= (64) 

We see that / plays, up to a factor, the role of 1/6 = GM^/Mp in the original model, and so 
the effective scalaron mass squared Mg is substituted here by a combination of squared masses of 
fermion and boson fields - confer Eq. (I61h . 
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With the anomaly-induced formulation we have recovered once more a “tempered” form of 
inflation, namely possessing an evolution pattern very similar to that indicated in Fig.[2l and thus 
presumably having a paved way to graceful exit. Not only so, now with the alternative formulation 
we have gained an additional bonus, to wit: the running Planck mass squared (I63p leads to a an 
effective gravitational coupling 


GN{a) 


Gn 

1 — / Ino ’ 


(65) 


where Gn = ^/Mp is the current value (a = 1). Let us assume that this situation is roughly 
maintained at low energies^. If we next insert (|65p into the general low-energy equation for energy 
conservation in the presence of variable cosmological parameters Gn and /?a, i.e. Eq. ([T]), we obtain 
a differential equation relating and pm- However, if we further assume matter conservation, 
then we can derive the evolution law for the vacuum energy density pA that is compatible with 
(|65P according to the Bianchi identity ([ 8 |). The exact result is a bit complicated, but in the limit 
I/I <C 1 (which is in fact the natural situation for a /3-function coefficient) it turns out that the 
final result for the matter-dominated epoch reads very simple |18] : 


p^{H) = pI + - Hi) , ( 66 ) 

where p\{Hq) = is the current value. This expression is noticeable as it shows that the vacuum 
energy of the anomaly-induced formulation of the Starobinsky inflation leads, at low energies, to 
a dynamical vacuum model of type Al, see Eq. (1141) . We can easily identify by comparison of the 
two equations that = //3. It shows that the coefficient v (which could have been introduced in 
(jl4|) on mere phenomenological grounds) can be interpreted as a /3-function coefficient for the CC 
running in semiclassical curved spacetime. One can indeed verify that the structure of / in (I 6 ip 
takes on the expected general form (|13p . This result reinforces the physical meaning of the general 
renormalization group equation of the cosmological term, Eq. (|10l) . in QET in curved spacetime. 

At the end of the day, we conclude that the class of vacuum models (jl 2 h encodes basic features 
of different realizations of Starobinsky’s type of inflation, mainly if realized in the anomaly-induced 
form. It remains to be seen what is the final observational status of the Starobinsky model. At 
the day of writing this work the situation is inconclusive [651178] and therefore we must wait for 
more input. 


7 Vacuum dynamics from ACDM in the current Universe 

In this section we go phenomenological and report on the confrontation of the dynamical vacuum 
models AI,A2, B1 and B2 defined in Sect. 3 against the latest observational data. All of these 
models can be solved both at the background and perturbative level. The Type-B ones, how¬ 
ever, are more difficult to deal with analytically and require more numerical work, specially in 
regard to the perturbation equations and the study of structure formation. Eor details see the 
recent works [221 ESI, which are to our knowlwdge the most comprehensive studies on these mat¬ 
ters currently existing in the literature. Here we limit ourselves to provide a quick summarized 
presentation of the solution of models Al and A 2 only, but we provide also some numerical results 
of the type-B ones. 

cannot exclude that as we run down to the low energy regime, corresponding to the present Universe, there 
might occur some additional, infrared, renormalization of the coupling /. This could enhance its value as compared 
to the UV regime. 
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Figure 3! The evolution of the energy densities (normalized with respect to Hq) for the type-A1 model d in the post¬ 
inflationary epoch until our days. The curves shown are: radiation (dashed line), non-relativistic matter (dotted line) and 
vacuum (solid line, in red). Used inputs: v = 10“^, = 0.27, = (1 + 0.227^1/) (A^^;,r2^,/i) ~ (3.04,2.47 x 

10“^/i.“^, 0.71). We have used SttG = 1 units. 


7.1 The cosmology of type-A models 

Let us consider the vacuum model A2 defined in Eq. (|16p . Using the equations of Sect. 2 we can 
solve for the energy densities as a function of the scale factor. For cold matter and radiation one 
finds; 

Pm{a) = pI, , pr{a) = (67) 

and for the vacuum energy density: 


pA{a) = Pa + Pm (C ^ - 1) (o - l) + Pr (C' ^ - 1) (a - l) , 

where we have dehned 




1 - 1 / 




1 - ly 


( 68 ) 


(69) 


1 - 3P/2 ’ ’ 1 - 2P 

Obviously, for P = 0 the model A2 becomes simply the basic running vacuum model A1 (introduced 
in [H]). The corresponding Hubble function reads 


H\a) = Hi 


1 + ^ - l) + ^ (a-^«' - 1 






(70) 


where the cosmological parameters satisfy the usual sum rule + 12)? + = 1. As expected, for 

ly = P = 0 (i.e. C = = 1) ths Hubble function H{a) takes on the form of the concordance model. 

In this case we recover pm ~ a~^, Pr ~ and, of course, also pA =const. As an illustration, 
in Fig. [3] the evolution of the matter and vacuum densities are plotted from the radiation epoch 
to our days. It corresponds to a type-Al model (P = 0) and we have assumed v = 10“^, i.e. a 
typical value obtained from the hts to the observational data, see below. We can see that the 
vacuum density performs a substantial evolution until it asymptotes to the current value pA —^ Pa- 
It certainly looks more natural than the situation in the ACDM model, where it stays p\ = p\ all 
the time since the radiation epoch. 
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Model V e x^/dof 


ACDM 

0.293 ±0.013 

- 

- 

567.8/586 

A1 

0.292 ±0.014 

±0.0013 ±0.0018 

- 

566.3/585 

A2 

0.290 ±0.014 

±0.0024 ± 0.0024 

- 

565.6/585 

B1 

u.za/ooi4 

- 

U.Uin o 013 

587.2/585 

B2 

U.OUU_Q gQ3 

—0 0039^^ ^^^^^ 

u.uuna_Q QQ 24 

u.uuna_o 0021 

583.1/585 


Table ll The fit values for the various models using SNIa+CMBH-BAO^z data, together with their statistical significance 
according to statistical test. For model A2 we provide the fit value of z/eff (see the text). For model B2, we have set u = e. 


To study the evolution of matter perturbations we generalize the standard second order dif¬ 
ferential equation m for the growth factor D = 5pm/pm for the case when the vacuum term is 
dynamical, pA 7 ^ 0. The result (after neglecting contributions which are subleading) is |22] : 

D + {2H + ^) D- (ATTGpm - D = 0, (71) 

with T = —pAlpm- For the ACDM model we have p\ =const. and hence T = 0, so that the above 
equation correctly reduces to the standard one [79] , 


7.2 Confronting the vacuum models A and B with observations 

Equation ([7T]) can be solved for type-A models (for which pm and pA have been given above) in 
terms of hyper geometric functions. In the case of type-B models the solution is more complicated 
and one has to use direct numerical methods |22l 123] , From the corresponding solution for D{a) 
we can e.g. study the linear growth rate of clustering [79], namely the logarithmic derivative of the 
linear growth factor D{a) with respect to In a: 


1 dD{a) d\iiD{a) 

D[a) din a dlna 


-(1 + ^) 


din D{z) 
dz 


(72) 


This quantity is important enough since it is measured observationally, so we can use it to in¬ 
vestigate the performance of our vacuum models. In Fig. [5] we compare the theoretical growth 
prediction derived from our fits to the data (see below) with the latest growth data (as collected 
e.g. by [80| and references therein). Specifically we plot the combined observable f{z)as{z), viz. 
the ordinary linear growth rate weighted by the rms mass fluctuation field, which it is claimed to 
have some advantages in practice [ 8 T]. The value of (Ts{z) can also be computed for each model 
through as{z) = as D{z)/D{0), where we use as = 0.829 ± 0.012 from Planck-l-WP [4]. 

It t is pretty clear from Fig. U] that the dynamical vacuum models under consideration are able 
to adjust the linear growth data in a way comparable to the ACDM. This was expected since the 
“running” of the vacuum density in our time does not depart too much from the constant p^ value. 
However, if one sets cq = 0 in the running models the departure from the structure formation data 
becomes dramatic. This is quite manifest in the anomalous behavior of the Cl curve in Fig.Uj 

We have performed the ht of the various vacuum models from the combined data on type la 
supernovae (SNIa) |82|, the shift parameter of the Cosmic Microwave Background (CNB) |1], and 
the data on the Baryonic Acoustic Oscillations (BAOs) |83|. In the last case we used the BAO^^ 
data based on the dz estimator defined in that reference - confer Table 3 in it. The basic htting 
procedure is explained in detail in Ref. |22] . Our results are collected in Table 1 above. Owing 
to some degeneracies in the two parameters of models A2 and B2, for the former we have set 
1 / = Vl2 (such that ^ = 1) and for the latter v = e. With this setting, and taking into account 
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f(z)crg(z) 




Figure 4) Comparison of the observed (solid points with vertical error bars) and theoretical evolution of the weighted growth 
rate f{z)crs{z) for the various A and B models. The uppermost (red) line corresponds to the ACDM model, used as a reference. 
The curve Cl (black line) that deviates significantly from the others in the two panels corresponds to model B2 for cq = 0. 
The curves have been obtained for the best fit values indicated in Tabled For more details, see m- 


that \u\, \v\ <C 1, we have from ([69]) ^ ~ 1 — t'eff) with t'gff = ^ The fitting value to this 

effective parameter is in fact the one quoted in Table 1 for model A2. 

Parameterizing the linear growth as f{z) ~ one can define the linear growth rate 

index m 7. It can be used to distinguish cosmological models, see e.g. [H]. For the ACDM 
model such index is approximated by 7 a ~ 6/11 ~ 0.545. For the dynamical vacuum models we 
make use of ()72p and the knowledge of the corresponding growth factor D{a) for each model, and 
we get 




(73) 


The growth rate index 7 for each model is defined from the value of the previous expression for 
z = 0. 

In Fig. [5] we plot the evolution of 7(2:) for the A and B type of vacuum models. In the same 
figures we can also see our determination of ^a{z) for the ACDM as a function of the redshift, and 
in particular we find 7a(0) — 0.58. From that figure we see that the growth index of the dynamical 
vacuum models (especially for type-A models) is well approximated by the ACDM constant value 
for z <1, while at large redshifts there are deviations. It is worth mentioning that the differences 
with respect to the ACDM are at the edge of the present experimental limits. For example, in a 
recent analysis it is found that 7 = 0.56 ± 0.05 and = 0.29 ± 0.01 [H]. The prediction of 7 for 
all our vacuum models lies within la of that range. Being the differences on the verge of being 
observed, we expect that in the future it should be possible to discriminate among the various 
ACDM models by this procedure. 

Finally, we briefly mention an alternative procedure to distinguish the running vacuum models 
from the standard ACDM. It is based on the cluster number counts method, see [6I1I60] for previous 
applications in this context. The method is based on the Press-Schechter formalism |85] and 
generalizations thereof [86] . We refer the reader to Refs. [221123] for a comprehensive presentation, 
and in what follows we show only the most important results. The main observable of the method 
is the fractional difference 5M{z)/M{z) in the number of counts of clusters between the vacuum 
models and the concordance ACDM model at any given redshift 2. The number of counts at 
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Figure 5: The evolution of the growth rate index, Eq. ( 1731 ) . for the various vacuum models under the same inputs 
as in the previous figure. 


redshift z is given by the following expression; 


Af{z) = 


47rr^(z) p{z) 
H{z) 


I f I da \ ^ , , ,,, 


(74) 


Here p{z) is the comoving background mass density. The original Press-Schechter function is [85] 
/psc(c) = -\/2/7r((ic/<7) exp(—(5^/2cr^). It depends on the parameter 5c-, the linearly extrapolated 
density threshold above which structures collapse. This parameter must be computed for each 
vacuum model using the non-linear perturbations equations [221 [23] . A generalized and improved 
form of the Press-Schechter function /psc(<7), which we adopt here, is the one provided by Reed et 
al. in [86]. Finally, r{z) in the above equation is the comoving radial distance out to redshift z, 
namely: 



(75) 


and a‘^{M, z) is the mass variance of the smoothed linear density field. It depends on the redshift 
z at which the halos are identified and is given by 


a‘^{M, z) = r k‘^P{k)W^{kR)dk . (76) 

Jq 

In this expression, P{k) is the power-spectrum and D{z) is the linear growth factor of pertur¬ 
bations. It is obtained from solving the perturbations Eq. ()7ip and is therefore characteristic of 
each model. Finally, W is the Fourier transform of the standard top-hat function with spherical 
symmetry - for details, see e.g. [22] , 

By selecting an appropriate redshift range for the observations, as well as a characteristic 
range of masses Mi < M < M2 for the observed clusters, one can provide definite predictions 
from (f74|l . Using the best fit values of the parameters of our vacuum models, as indicated in 
Table 1, and the full machinery of the generalized Press-Schechter formalism, we can compute the 
redshift distribution of clusters predicted in each vacuum model and subtract the corresponding 
prediction from the ACDM, i.e. the quantity 5M(z) = M{z) — M{z)\cdm ■ Then we can compute 
the fractional difference 6M(z)/M(z) \cdm , which encodes the relative deviations with respect to 
the concordance model. The numerical results are shown in Fig.[6l 

The vacuum models are seen to be clearly separated by the cluster number counts method, 
with the type-B ones providing an excess in the number of counts, and the type-A models a 
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6N/M, ( 10''> M h/Mg >10’="^) 



Figure 6: Comparison of the fractional difference 5M/M in the redshift distribution of cluster number counts (in the 
indicated range of masses) for the vacuum models under consideration with respect to the concordance ACDM model, i.e. 
&M = M — AfACDM- Inputs as in the previous two figures. 


defect, as compared to the ACDM. The deviations can be significant (+30 — 40%) for the former 
(especially for B2) but are moderate (say —10% to —20%) for the latter, at around the optimal 
redshif z ~ 1.5 where it lies (approximately) the maximum total number of counts (cf. [22]). The 
fractional differences are still bigger for higher redshifts, reaching +50% for B2 at 2 ; ~ 2. Beyond 
redshift z = 3, however, the statistics of the total number of counts depletes significantly and as a 
result this range is not appropriate for practical measurements. 

8 Conclusions 

In this work we have discussed various aspects of the “ running ACDM cosmology” (a large class of 
dynamical vacuum models) in which the vacuum energy density can be expressed as a power series 
of the Hubble function and its cosmic time derivatives. Some of these models can be well motivated 
within the context of quantum field theory (QFT) in curved spacetime and they can provide an 
overarching description of the cosmic history, namely starting from the early inflationary times till 
our dark energy (DE) days. For the study of the current Universe the series naturally terminates 
at the level of the and H terms, but the higher order ones are crucially important for a proper 
description of the inflationary phase. 

The class of ACDM models involving higher powers of the Hubble rate (typically have 
a non-singular starting point characterized by an initial de Sitter period of rapid inflation, pre¬ 
sumably triggered by the high energy dynamics of a Grand Unified Theory at a scale Mx near 
(but below) the Planck scale. Explicit solution of the model shows that the inflationary period is 
followed by successful “graceful exit” into the standard radiation epoch. The Universe continues 
into the cold dark matter era and finally leads to a dark energy epoch with a very small amount 
of vacuum energy. 

It is important to emphasize that the post-inflationary period leads to a cosmological regime 


27 













that is very similar to the concordance ACDM model, but with a distinctive feature: the vacuum 
energy is not rigid but mildly dynamical ~ cq + u H"^). This feature can be tested and may 
prove a low-energy “smoking gun” of the underlying vacuum dynamics. Let us remark the recent 
indications on dynamical DE in the current data m- The slow vacuum evolution in our recent 
past, still going on in our days, could be responsible for the dynamical DE and could provide an 
alternative explanation to it which is completely different [22] from quintessence approaches and 
the like. Similarly, the same unified vacuum structure explaining the DE is responsible, at the 
very early times, of the inflationary period and without invoking at all any sort of ad hoc inflaton 
fields. Therefore, in contrast to the ACDM model, the very early times of the cosmic history in 
the ACDM cosmology is characterized by a powerful dynamical vacuum (~ ~ ^x)i which 

triggers inflation and upon decay into matter becomes much more moderate ( H^> Hq) near 

our present until effectively behaving as DE. 

Quite noticeably, the running ACDM model provides a natural explanation for the huge entropy 
of the current Universe. It also emerges from the primeval vacuum decay in the early Universe. 
While the details of the vacuum decay depend on the Grand Unified Theory that brings about 
inflation, the final prediction of the entropy is universal and does not depend on the particular 
GUT implementation. The mechanism we have described here explains the value of the entropy in 
the current horizon. No horizon problem exists at all in the ACDM model since all of the points 
of the current Hubble sphere remain causally connected as of the early times when a huge amount 
of relativistic particles emerged out of the decaying dynamics of the primeval vacuum. 

We have signaled some potential connection of the ~ structure of the AGDM cosmology with 
Starobinsky inflation, characterized by B? ~ H^. The two sorts of models present some similarities 
but they actually involve also important differences, which depend on particular implementations. 

After the inflationary epoch is left behind, the total amount of vacuum energy decreases very 
fast, the primordial nucleosynthesis can operate normally within a standard radiation epoch and 
the Universe can go through the cold dark matter era until the present DE era. The latter appears 
here as a slowly time-varying vacuum-dominated epoch. The low-energy vacuum dynamics can, 
however, adopt different forms (which we have called type-A and type-B). 

We have shown that at leading order all these forms are admissible structures for the consistent 
description of our Universe, but they exhibit some differences that can be checked observationally. 
In particular we have confronted our vacuum models against the Hubble expansion data and struc¬ 
ture formation. In addition we have assessed their considerably different capability in populating 
the Universe with virialized structures at different redshifts as compared to the ACDM model. 

The current Universe appears in all these models as FLRW-like, except that the vacuum energy 
is not a rigid quantity but a mildly evolving one. The typical values we have obtained for the 
coefficients u, a and e responsible for the time evolution of pA lie in the ballpark of rsj 10 " 3 . This 
order of magnitude value is roughly consistent with the theoretical expectations, some of them 
interpreted in QFT as one-loop /3-functions of the running cosmological constant [UliHlIIg]. It 
is rather encouraging since it points to a fundamental origin of the theoretical structure of the 
ACDM models in the context of QFT in curved spacetime. 

To summarize, the class of ACDM cosmologies may offer an appealing and phenomenologically 
consistent perspective for describing inflation and dynamical dark energy without introducing 
extraneous dark energy fields. Ultimately they might offer a clue to better understand the origin 
of the A-term and the cosmological constant problem in the context of fundamental physics. It 
would be a timely achievement, if we take into consideration that we are currently approaching 
the centenary of the introduction of the cosmological term in Einstein’s equations in 1917 [T]. 
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